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1. Introduction 

In the last decade, percolation has enjoyed the attention of conformal field theorists who 
have sought to calculate crossing probabilities for various aspect ratios and geometries. 
In rough terms, calculations generally involve mapping percolation to a 1-state Potts 
model, constructing a correlation function corresponding to the boundary conditions 
necessitated by the crossing cluster of interest, and finding a differential equation thereof. 
In 1991, Langlands et al [2] were the first to investigate "the universality of crossing 
probabilities in two-dimensional percolation". Shortly afterwards, Cardy [3] obtained 
an exact equation for the probability of a crossing cluster on a rectangle for different 
aspect ratios, using conformal field theory. 

In 1996, Watts extended Cardy's results to obtain an exact equation for the 
probability of a cluster crossing both horizontally and vertically [4]. In the same 
year, after numerical work by Hu and Lin [5], Aizenman proved that the probability 
of more than one crossing cluster is finite in the thermodynamic limit [6]. This does 
not contradict the rigorous result that the number of infinite clusters is either 0, 1 or 
oo with probability one [7, 8]. The asymptote for the probability to obtain n distinct, 
simultaneously crossing clusters has been calculated by Cardy [9], who also considers 
spanning on a cylinder. 

We define all these probabilities systematically below. In this paper, we provide 
accurate numerical data for the many analytical results now available, many of which 
are verified for the first time. We also provide data, particularly for wrapping on the 
cylinder, for which there is currently no theoretical explanation. After briefiy discussing 
the results for some exotic cluster configurations, we present and discuss results for 
crossing, spanning, and wrapping probabilities on rectangles and cylinders of various 
aspect ratios for site and bond percolation. The results are of such accuracy that we 
are able to make firm statements as to the validity of various estimates, conjectures 
and formulae made in the past. An appendix collects together some technical notes on 
identifying various cluster types when simulating percolation. 

2. Observables 

In this section we give the various observables that we measured for site and bond 
percolation. We consider a square lattice in which sites are linked via bonds. To avoid 
confusion we stress that boundaries are made of sites only [2]. In site percolation sites 
are occupied with probability p^''' and all bonds are active. In bond percolation all sites 
are occupied, while bonds are active with probability p^^\ Two sites belong to the same 
cluster if they are both occupied and if there is a path between them along active bonds 
and occupied sites. A cluster is the set of sites belonging to it. 
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Figure 1. A realisation of a lattice which is consistent with figure 8. As indicated, 
the clusters are of type S^jv,£;}i S|5 j^/j and Sj^vi'}. The hatched square in the centre 
will be filled to give rise to the border configuration shown in figure 9. 



2.1. Open boundary conditions 

Each cluster can be characterised by the borders it touches. If the borders of the lattice 
are labeled iV, E, S, W, as in figure 1, then each cluster is assigned a subset of these 
labels, indicating the borders it touches. There are 15 different combinations of cluster 
labels in which at least one border is touched. In the following, these combinations are 
called "types", Sq. A cluster of type Sq touches all borders in the set fl. For example, 
the three clusters shown in figure 1 are of type S^at^^;}, Sj^^iy} and S^w}- 

We distinguish between clusters touching only a set of borders and at least a set of 
borders. A cluster of type Sq is also of type Uq' if ^2 ^ fi'. For example, a cluster of 
type Sjs^iy} is also of type ^{s}, ^{w} and ^{s,w}- A "crossing cluster", i.e. a cluster 
that connects two opposite borders, is therefore either of type UjAr^g} or U{£;,vk}4 

Normalised histograms 'PAr(T, n, r) have been generated for all 15 S and U types, 
where is the size of the lattice (number of sites) and r is its aspect ratio (length over 
height). The histogram estimates the probability that a random realisation contains n 
clusters of type T. It is worthwhile pointing out that even though S and U histograms 
are correlated, they cannot be derived from each other, because the number of clusters 
of a particular U-type must be determined on a per-realisation basis. 

The moments associated with each histogram are defined by 

oo 

MN{T,m,r) = Y.VN{T,n,r)n"' . (l) 

n=0 

A hat on a quantity, for example Ai, indicates its expectation value, i.e. the value in 
the limit of an infinitely large ensemble and system size. 

2.2. Cylinder 

In addition to the type classification above, we have identified wrapping and spanning 
clusters on a cylinder. Our convention for cylinders is illustrated in figure 2, which 

if In the following, we will apply the term "crossing" only to systems with open boundaries, and reserve 
the terms "spanning" and "wrapping" for cylindrical boundary conditions. 
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Figure 2. A realisation of percolation on a cylinder that would contribute to 
7'Ar(0+, 1, r) and PAr(l+, 1, r), but not to VNiO, 1, r) or Pjv(l, 1, r). 

defines the aspect ratio as circumference over height. A wrapping cluster is then a 
cluster which winds around the cylinder, i.e. it provides a path of length of 2tt. A 
spanning cluster, meanwhile, connects the bottom and top of the cylinder. 

The resulting histograms are PAr(l, n, r) for the number distribution of clusters that 
only span (rather than wrap) the cylinder, and 'PAr(0,n, r) for the number distribution 
of clusters that only wrap (rather than span) the cylinder. 

In addition, Vn{\^i 1,?") counts the number of realisations with a single spanning 
cluster that may also be wrapping. Similarly, VNiO^Ai"^) counts the number of 
realisations with a single wrapping cluster that may also be spanning. For an example, 
see figure 2. The distinction between I and 1^ (0 and 0^) disappears for more than 
one spanning (wrapping) cluster, because the existence of two or more simultaneously 
spanning (wrapping) clusters prohibits wrapping (spanning) at the same time. 

3. Results 

We have generated the histograms mentioned above for site (p^°' = 0.59274621 [10]) 
and bond percolation [p^^^ = 1/2 [11]), each for three different system sizes, = 
30000^, 3000^, 300^. If not mentioned explicitly, the results presented are for A^ = 30000^. 
Henceforth, site percolation is indicated by a superscript (s), and bond percolation by 
(b). In each of these simulations, 14 different aspect ratios were realised while keeping 
the area constant [12]. These aspect ratios are: 30/30, 36/25, 45/20, 50/18, 60/15, 
75/12, 90/10, 100/9, 150/6, 180/5, 225/4, 300/3, 450/2 and 900/1. Three different 
boundary conditions were applied to each aspect ratio, corresponding to an open system, 
a cylinder glued vertically, and a cylinder glued horizontally. "Gluing" is our technical 
term here referring to the procedure for applying periodic boundaries. For the definition 
and explanation of these technicalities, see [1] and the appendix. The random number 
generator used is described in [13]. Between 10^ and 2 x 10^ independent samples 
were produced for each system size and percolation type, using up to 61 undergraduate 
computers when idle. 
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The different aspect ratios are derived from the same set of patches (s. appendix), 
by gluing them together to form rectangles in random permutations and orientations 
(for details see [1]), and are therefore not statistically independent. However, their 
correlations are assumed to be very small. When considering the outcome for several 
aspect ratios at the same time, one could multiply the error bars by the square root 
of the number of aspect ratios considered, as if each aspect ratio were based only on a 
subset of the original sample. Because this procedure assumes maximum correlations 
between the patches, even though they are randomly permuted, rotated and mirrored 
between different aspect ratios, this is a strong overestimation of the correlations. 

We group our results as follows: 

• Corner clusters, types S{Ar,£;}, S{£;,5}, S{5,h/} and S{n,w}- 

• "Three-legged" clusters, types S{n,e,s}, ^{e,s,w}, ^{n,s,w} and S{n,e,w} 

• Types U{Ar^5} and V^[E.yv}^ the subject of Cardy's predictions [3, 9], as well as 
spanning and wrapping clusters on a cylinder. 

• Type S{N,E,s,w}^ in answer to Watts' prediction. 

Apart from some straightforward arguments, there are no theoretical predictions for 
corner and three-legged clusters, which we will refer to as "exotic" types. 

3.1. Exotic types 

3.1.1. Corner Clusters. For completeness we include the following results for corner 
clusters, i.e. types SjAr^^}, S{£;^5}, SjAr^iy} and S{5^vi/}- In a rectangle, all corners 
are equivalent. Naively one expects these clusters to be arranged like onion skins. 
Assuming scale invariance, this suggests a logarithmic dependence of their average 
number Al7v(S, l,r) on the lattice size. This, however, is not the case, as can be seen 
in figure 3, which should show collapsing lines if N enters only as a factor In(A^). It is 
clear that the graph must level off for r ^ oo, but this region is not yet reached, even 
for r = 900 and any A^. The lines are remarkably straight, but they seemingly cannot 
define a universal exponent. 

In the Ising model, the corner magnetisation has been calculated analytically by 
Davies and Peschel [14] using a corner transfer matrix approach [15]. A similar approach 
seems to be suitable for corner clusters in percolation. 

3.1.2. Three-legged clusters. Three-legged clusters touch three borders in a T-like 
manner. There are four different types of three-legged clusters, but the symmetry of a 
rectangle splits them into two pairs, while a rotation of the system by 7r/2 transforms r 
into 1/r so that 
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Figure 3. MN{S{M,E}A,'r) / \n{N) for N = 3002,3000^,300002 and site and bond 
percolation. 



The statistics of these clusters turns out to be universal. This is possibly not very 
surprising, since they involve a crossing path (see below). However, what is surprising is 
the asymptotic value of 'P(S{7v £; ^j, l,r) in the large r limit, when the vertical crossing 
probability approaches 1. At r = 100 and = 30000^, these asymptotic values are 
0.5004(4) and 0.5002(3), for site and bond percolation respectively. At r = 900, the 
values are 0.4999(4) and 0.5004(3). Analysing smaller system sizes reveals that this 
quantity is very sensitive to finite-size effects: while even the smallest system size 
reaches a value very close to 0.5 at r ^ 6, for small A^ deviations towards higher 
probabilities occur at large r. Evidently, the probability of a three-legged cluster tends 
to the occupation probability of sites on the lattice for r ^ N, where the system 
effectively becomes a one-dimensional strip of length A^ and height 1. The region where 
the probability remains close to 0.5 becomes larger as the system size increases, and for 
A^ = 30000^ we were unable to detect any significant deviation from 0.5 for r > 9. We 
therefore conjecture that 

}^nS{N,E,s}A,r) = 1/2 . (2) 

This can be understood as the probability of intersecting a vertically crossing cluster 
when vertically cutting a long, narrow percolating system. However it remains unclear 
how to derive this limit analytically. 
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3.2. Results related to Conformal Field Theory 

We now present numerical results that are related to theoretical predictions based 
on conformal field theory. This includes measurements of many formerly unknown 
quantities, conjectures, and comparisons with exact results, which all give further 
support to the conformal invariance of critical percolation. 



3.2.1. Crossing probability with open boundaries. Cardy's seminal paper [3] contained 
a comparison between his exact result for the crossing probability and the numerical 
results obtained by Langland et al [12]. These were based on systems with = 200^ 
sites and r = 1---7.35. Later studies by Shchur and Kosyakov [16, 17] investigated 
"PatCJiat^s}, n, r) at r = 1 for n > 1 with very small systems, < 64^. Shchur [18] 
later extended these results to systems up to size N = 256 x 3200, apparently still 
encountering finite-size corrections. Other studies, such as by Sen [19, 20] and Hove 
and Aharony [21], used similar system sizes, while also considering other properties of 
spanning clusters. 

Using the data presented in this article, it is possible to compare Cardy's prediction 
with much greater accuracy, based on systems with A^ = 30000^ sites and r = 1 ■ ■ • 900. 
However, for very large (and very small) values of r, relevant clusters become either too 
rare to give any reasonable estimate for the associated probabilities, or their number 
count becomes too broadly distributed. 

The crossing probability is the probability to find at least one crossing cluster in a 
particular direction. By symmetry 

V{\J{N,s}, n, r) = V{\J{E,w}, n, r"^) (3) 

which has been used in the data presented below. Consequently, results for r > 1 and 
r < 1 are not statistically independent, since they are based on the same realisation 
(but supposedly different clusters contribute). Using the short-hand notation 

oo 

V{J,>n,r)=Y,V{J,m,r) (4) 

m=n 

for arbitrary cluster type T, Cardy's exact result reads 

V{U{^,s}, > 1' = ^ (^' I ^' ^) 

where 

K[u) is the complete elliptic integral of the first kind and 2-^1 is the hyper geometric 
function. 

Figure 4 shows the difference between the numerical result and the exact value 
from (5) in units of standard deviations. From this plot it is clear that the systematic 
deviation for large r observed in [3, 12] was only a finite-size problem. 
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Figure 4. Deviation between analytic (5) and numerical results for the reduced 
crossing probability 'PAr(U{jv.s}i ^ li in units of standard deviations of the numerical 
results versus aspect ratio r, for bond (circles) and site (triangles) percolation with 
N = 30000^ sites. 



The asymptotic number distribution of crossing clusters in percolation with open 
boundaries has been derived analytically [9] 

V{\J{N,s},n,r) ^ Cx{n)exp(^~7in(^n-^^ /r^ for n > 1 (7) 

in the limit of ^ oo and r ^ 0. In general, the amplitude Cx{n) is not known exactly, 
but it is universal and can be derived from numerics by fitting (7) against the numerical 
data in an appropriate region of r values for each n separately. The range of aspect 
ratios used in the fit is determined by two competing interests: The fit should include 
as many points as possible, but exclude aspect ratios where the asymptotic behaviour 
has not yet set in. This determines the largest r in the fit. The smallest r is given by 
the value of r for which n simultaneously crossing clusters are observed at least once. 
Fitting ranges are given in all tables below. 

In contrast to spanning on the cylinder, Cardy's asymptotic formula (7) does not 
distinguish between a cluster crossing in exclusively one direction and a cluster crossing 
in possibly more than one direction. For consistency in our notation, Cx{n) refers to 
amplitudes of exclusively vertically crossing clusters. However, at n = 1 (7) applies 
to vertically crossing clusters irrespective of other clusters, to which we assign the 
amplitude C^{1). For n = 1 Cardy's prediction (7) therfore reads 

P(U{^,s},l,r)-C+fl)exp(-^7r/r) . (8) 
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There is no prediction for exclusive crossing, ^(UiAr^^}, 1, r) — P{}^{n,e,s,w}A^'^)i 
however, it is consistent to fit the latter against 

Cx(l)exp(-a(f)/r) . (9) 

The fits throughout the article depend to some degree on the choice of the fitting 
interval. A priori, it is unknown where the asymptotic behaviour sets in (in the sense 
that the deviation from the asymptote is smaller than the numerical error). It is not 
possible to determine whether a deviation of the numerical results from the analytical 
value is due to statistical fluctuations or due to a wrongly chosen fitting interval. The 
errorbars given can only reflect the former. However, the error indicated should include 
the exact result if it is fitted against the corresponding function in the range given. 

Figure 5 shows the numerical data for the probability of finding n = 1,2,3,4 and 
n > 1 vertically crossing clusters in reduced form, i.e. hi{V / {1 — V)). The data are fitted 
to the asymptotic formulae (7), (8) and (9). As discussed below, the latter coincides 
asymptotically with (12), which is shown as a dotted line. For completeness, data for 
7^(U{7v,s}5 ^ l;''^) are shown together with the exact result (5). 
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Figure 5. Numerical data {N = 30000^, bond percolation, open boundaries) for 
the reduced probability ln('P/(l — V)) oi n = 1,2,3,4 and n > 1 vertically spanning 
clusters versus aspect ratio r. The data n > 1 refer to 7^(U{Ar^5},> l,r). The long 
dashed lines are fits according to (7) ((8) and (9) for ti = 1) and Tab. 1, while the full 
line gives the exact result (5). The dotted line gives the exact result (12). 



It is clear that, asymptotically, the relative difference between the probability of 
more than one and exactly one crossing cluster, 
P(Uws},>l,r)-P(U{^,5},l,r) 
nU{jv,5},>l,r) + ^(U{^,5},l,r) 
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vanishes, and the crossing probabihty becomes dominated by the probabihty of a single 
crossing cluster. Therefore, the amplitude (1) in (8) is known exactly, namely 

C5(l) = 2^^' ^ 1/3)2 = 1-4263482556253 • ■ • (11) 

using the expansion of (5) provided by Ziff [22, 23]. 

Equation (10) implies that the probability of exactly one crossing cluster exclusively 
in one direction, 'P(U{Ar^5}, 1, r) — V{\}{n,e,s,w}A^'^)^ is the dominating term in the 
difference 

^(U{JV,5}, > 1, r) - V{{}{N,E,sm^ 1' 
rj / 4 5 



r(i/3)r(2/3) ^^H^'^'s'^'s'^'j ' ^^^^ 

with rj as in (6). This relation is derived from (5) and (29), and for small r can be 
expanded as 

r(l/3)r(2/3) ^^P(~^/^)(^ - ^exp(-7r/r) + ^exp(-27r/r)---) (13) 

which is based on Eqn. (16) in [22] and the definition of the generalised hypergeometric 
function. Therefore the two parameters appearing in (9) are known exactly: 

16 

Cx(l) = r(i/3)r(2/3) = 4.41063116337433639 ■ • • (14) 

and 

a{l) = n (15) 

which corresponds to n = —1 or n = 3/2 in (7). 

The numerics are in very good agreement with the exact results. The estimate 
for Cx (1) shown in Tab. 1 agrees perfectly with (11), while Cx(l) has a surprisingly 
large error. However, the result still covers the exact value (14). The error is due to a 
narrow fitting range, forced by the late onset of asymptotic behaviour. A larger fitting 
range gives a much smaller error and reduced goodness-of-fit [24]. For a(l) it is found 
numerically that a^^^l) = 3.18(6) and a'°^(l) = 3.15(5), in perfect agreement with (15). 

The amplitudes listed in Tab. 1 have also been used to plot the dashed lines in 
figure 5. They fit an exponential fairly well (using the l"*" value in Tab. 1): 

C^\n) ^ 0.5422(10) exp(0.9619(14)n) (16) 
C^=^(n) ^ 0.5434(9) exp(0.9617(12)n) . (17) 



3.2.2. Spanning probability on a cylinder. Cardy has also given asymptotes for 
spanning events on a cylinder [9]. These clusters have been investigated numerically 
several times [21, 16, 18]. However, this work uses system sizes three orders of magnitude 
larger than in former studies. The probability of obtaining n distinct spanning clusters 
on a cylinder is expected to behave in the limit of small r and large L like 

P{\,n,r)^Q{n)exp{-'^7T{n'-^)/r) . (18) 
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Table 1. Crossing with open boundaries. Amplitudes are defined as in (7-9), derived 
from numerical simulations of bond {C^^) and site (C^ ') percolation with N = 30000^ 
sites. The r range defines the fitting region (see section 3.2.1). The + marks , i.e. 
not an exclusively single crossing cluster in one direction. 



n 


r range 


bond 




r rang 


e site 




1+ 


9/100- 


-25/36 


1.423(2) 


9/100 


---25/36 


1.426(2) 


1 


9/100- 


-20/45 


4.811(629) 


9/100 


---25/36 


4.560(510) 


2 


20/45- 


-36/25 


3.530(12) 


20/45 


---36/25 


3.553(10) 


3 


30/30- 


-45/20 


9.608(37) 


30/30 


---45/20 


9.599(32) 


4 


45/20- 


-60/15 


27.641(161) 


45/20 


---60/15 


27.658(140) 



Table 2. Spanning on a cylinder. Amplitudes are defined as in (18), where + marks 
the value of C+(l) (see (19)). 



n 


r range 


bond 




r rang 


e site 




1+ 


10/90- 


-30/30 


1.2217(4) 


10/90 


---30/30 


1.2222(4) 


1 


9/100- 


-30/30 


1.7198(11) 


9/100 


---30/30 


1.7225(10) 


2 


25/36- 


-36/25 


5.1829(256) 


25/36 


■--36/25 


5.2105(218) 


3 


36/25- 


-50/18 


15.1212(764) 


30/30 


---50/18 


15.0227(649) 


4 


45/20- 


-60/15 


45.0059(3280) 


45/20 


---60/15 


44.5445(2780) 



The existence of a wrapping cluster prevents more than one spanning cluster, see 
figure 2. For n = 1 one can distinguish between exclusively spanning clusters and 
spanning clusters that may also wrap. Allowing for wrapping clusters Cardy predicts 
the asymptote 

P(lM,r)^C+(l)exp(-^7r/r) . (19) 

The numerical results are in full agreement with (18) and (19). The corresponding 
amplitudes are shown in Tab. 2. The fact that the amplitude in (19), C'^{1), is 
slightly smaller than the corresponding amplitude in (18), C|(n), does not contradict 
P(l, l,r) < V{\^^ since the exponentials differ by a factor exp(7r(7/24)/r) in favour 
of spanning without restrictions on wrapping. It would require an r > 2.67 to suppress 
this factor enough to equalise both probabilities. However, at such large values of r, 
Eqn.'s (18) and (19) are not valid any longer. 

The amplitudes fit an exponential extremely well (where for n = 1 wrapping was 
not allowed, i.e. the result l"*" in Tab. 2 was ignored): 

Cl^^n) ^ 0.5788(11) exp(1.0892(17)n) (20) 
C(=)(n) ^ 0.5816(10) exp(1.0860(14)n) . (21) 

The resulting plot of the reduced probabilities looks very similar to figure 5, so we 
omit it here. 
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Table 3. Wrapping on a cylinder. Results for the fits of 7^(0, n,r) to (23). The r 
range defines the fitting region, where the right-hand value is for the largest aspect 
ratio in which n wrapping clusters occured at least once. The + marks the value of 
Cq (1) that allows a spanning cluster. 



n 


r range 


(bond) 




ag>(n) 


r range 


(site) 






1+ 


25/36 • • 


• 100/9 


1.0654(12) 


-1.0755(10) 


25/36 • • 


•100/9 


1.0634(10) 


-1.0739(9) 


1 


25/36 • ■ 


•60/15 


1.7734(82) 


-3.0536(52) 


25/36 • • 


•60/15 


1.7712(70) 


-3.0522(44) 


2 


15/60 • ■ 


•36/25 


1.8134(55) 


-6.9320(92) 


15/60 • ■ 


•45/20 


1.7995(46) 


-6.9118(78) 


3 


12/75 • ■ 


•25/36 


3.6128(216) 


-17.8021(326) 


12/75 • ■ 


•30/30 


3.6600(186) 


-17.8764(278) 


4 


10/90 • ■ 


• 20/45 


6.6734(609) 


-32.6919(752) 


10/90 • • 


• 20/45 


6.6533(520) 


-32.7071(645) 



3.2.3. Wrapping on a cylinder. Aizenman's original statement [6] regarding the 
number of crossing clusters can also be applied to the number of wrapping clusters 
on a cylinder. In the limit of large r 

ln(P(0,n,r)) G C(rn2) (22) 

according to a hand- waving scaling argument by Cardy [9] . There is no better estimate 
known, so by fitting each set of histograms to 

p(0, n, r) = Coin) exp(ao(^)r) (23) 

one can determine the n dependence of ao(n) and the amplitude Coin). Tab. 3 
shows the corresponding results. It turns out that ao(^) fits very well a second order 
polynomial with coefficients 

"o^H ~ -3.150(11)^2 + 5.51(4)^-5.41(3) (24) 
"oH ~ -3.176(10)n2 + 5.61(4)r2 - 5.48(3) . (25) 

The amplitude Coin) again fits an exponential, but only if the very first value, Co(l), 
is neglected. In this case we find 

C!o\n) ^ 0.483(5) exp(0.662(4)n) (26) 
eg' (n) ^ 0.473(4) exp(0.670(3)n) . (27) 

We stress again that the ambiguity in the choice of the fitting ranges introduces an 
error which is not refiected in the numerical error given. We find it therefore justified 
to conjecture that 

7 

ao{n) = -nn^ + -7r(ra - 1) (28) 

which is somewhat surprising, since we naively expected similar arguments [9] to those 
for (18) should apply, giving rise to a leading term — Ivrn^. The numerical results are 
shown together with the proposed analytical behaviour according to (23) using the data 
in Tab. 3 in figure 6. 
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Figure 6. Numerical data {N — 30000^, bond percolation) for the reduced probability 
ln(P/ (1 — V)) of n = 1, 2, 3, 4 and n > 1 wrapping clusters on a cylinder versus aspect 
ratio r. The data n > 1 refer to V{0, > l,r) and look qualitatively similar to those 
shown in figure 5. The long dashed lines are fits according to (23) and Tab. 3. 



3.2.4- Spanning simultaneously in both directions. Watts [4] has exactly calculated the 
probability of a cluster that crosses both directions simultaneously in a system with 
open boundaries, V{\J{n,e,s,w}, l:^)^ given by 

V{\J{N,E,S,W}^^r)= (29) 

TxlF^^ H3'3'3'^j~f(i)rx|) ' 

with rj as defined in (6) and 3-F2 being the generalised hypergeometric function. The 
first term on the right-hand side of (29) is identical to 'P{\J{n,s},^ ^,^)- According 
to the expansion provided by Ziff [22] and consistent with (8), for small r this term 
is proportional to exp( — |7r/r), while, according to (13), the second term decays even 
faster, namely exp(— vr/r). This is of course what one expects, because the probability 
of spanning in both directions is, away from r = 1, dominated by the probability of 
spanning in the longer direction. 

Figure 7 shows the difference between our numerical results and Watts' prediction. 
As in figure 4, the plot suggests that the deviations observed in [4] are only finite-size 
corrections. 



Crossing, spanning and wrapping in two-dimensional percolation 



14 




Figure 7. Deviation between the analytic result (29) and the numerical results for 
the reduced crossing probability 'Pn{^{n.e.s,w}i ^^f)! in units of standard deviations 
of the numerical results versus aspect ratio r, for bond (circles) and site (triangles) 
percolation with N — 30000^ sites. 



3.3. Finite-size corrections 

Since the system sizes investigated are huge compared to former studies [12, 16, 18], one 
might be inchned to completely ignore finite-size corrections. In order to estimate their 
strength. Tab. 4 lists numerical results of various quantities for different system sizes and 
compares them with estimates for the value of these quantities in the thermodynamic 
limit found in the literature or their exact results. For large systems, our numerical 
results agree very well. For smaller systems, there may be some mild corrections. From 
the three-legged cluster data (sec. 3.1.2) we expect that finite-size corrections may 
become visible if one side of the rectangle has a length < 300, which is not the case for 
any aspect ratio we have simulated in a system with = 30000^ sites. 

The results also indicate that the estimate of p^"' by Newman and Ziff [10] is 
valid within numerical error. This is corroborated by figure 4 and figure 7, where 
site percolation does not seem to show stronger deviations than bond, for which p^^^ is 
known exactly. 

4. Conclusion and Discussion 

Our numerical results represent probably the most comprehensive and most up-to-date 
study of crossing, spanning and wrapping probabilities. We have presented results for 
"exotic" cluster types, for which an analytical description is still lacking. With regards 
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Table 4. Comparison of our numerical values and estimates found in the literature, 
for different system sizes. "Estimates" are values for those quantities cited in the 
literature in the thermodynamic limit. 



Quantity 

7^Ar=30002(U{Ar,5}, > 1,1) 
7^Ar=300002(U{iv,S}, > 1, 1) 

PjV=30002(U{jv,S}, > 2,1) 
■pAr=300002(U{Ar^5i,, > 2, 1) 

'PAr=30002(U{jv,_E,S,lV}) I7 1) 
'PM=3mm^{^{N,E,S.W}i 1) 1) 

:P]V=3000^(I,> 1,1) 
7^^=300002(1, > 1, 1) 



site bond 



0.4995(4) 0.5012(4) 

0.4998(3) 0.4999(4) 

6.58(6) X 10^3 6.69(7) x 10"'^ 

6.71(6) X 10-3 6.66(6) x lO^^ 

0.3221(4) 0.3230(4) 

0.3219(3) 0.3226(4) 

0.6360(4) 0.6368(4) 

0.6361(3) 0.6364(4) 



Exact value /Estimate 



1/2 exact [2] 



6.58(3) X 10-3 [16] 



0.322120455- ■• exact[4] 



0.63665(8) [21] 



to the predictions of conformal field theory, our numerical data affords a comparison 
between numerical and analytical results, and we give further support to conformal 
invariance in percolation and other critical phenomena [25]. We believe that deviations 
from the predicted behaviour observed in the literature are most likely due to finite-size 
effects. 

We have also calculated the amplitudes listed in Tab. 1, Tab. 2 and Tab. 3, 
which might be of help to theorists. For three-legged clusters, we have conjectured 
the asymptote (2), for wrapping probabilities the form (28). 
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Appendix 

Here we list a few technicalities on how to identify cluster types, etc. The identification 
methods are applicable in any simulation of percolation where the boundary can be 
represented in the Hoshen-Kopelman [26] form. 

Simulation Method 

Our method performs asynchronously parallelised percolation on distributed machines 
[1]. In principle, the method relaxes all the standard constraints in numerical simulations 
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Figure 8. An example of a labeled border. A negative entry indicates a cluster size, 
a positive entry is a pointer to another site, and a zero entry indicates that the site 
is not occupied. The configuration shown is a border representation of a full lattice 
as in figure 1. The dashed pointer applies for periodic boundary conditions along the 
vertical borders, joining the two highlighted sites. 

of percolation, such as CPU power, memory, and network capacity. It is especially suited 
for calculating cluster size distributions and finite-size corrections, crossing probabilities, 
and, by applying the corresponding boundary conditions, distributions of wrapping and 
spanning clusters on different topologies, e.g. cylinder, torus, or the Mobius strip. 

The method is based on a master /slave parallelisation, where slaves send "patches" 
(specially prepared borders representing the lattice) to a master node, which "glues" 
these patches together. After a path compression, which is essentially a form of 
"Nakanishi label recycling" [27, 28] where bulk sites are considered inactive, the result 
is a single border as shown in figure 8. The term "gluing" will be used to indicate 
that a configuration is updated to account for a link introduced between two sites 
or two boundaries. The Hoshen-Kopelman (HK) algorithm [26] provides the data 
representation, which is the key of this method. Extremely large system sizes can 
be simulated. 

An example of a border representation is shown in figure 8 (in the following all 
examples are based on site percolation). Each site on the border is indexed in a clockwise 
manner starting with 1 in the upper left corner. Each site also contains a label. If the 
label is the site is not occupied. If it is positive it is interpreted as a "pointer" to 
the index of another site in the same cluster. If it is negative, it is called a "root" , and 
the magnitude of the negative number indicates the size of the cluster the site belongs 
to. The HK algorithm ensures that all sites of a cluster form a single tree, with a root 
that carries a negative label. In this way it is possible to identify the cluster each site 
belongs to by its corresponding root. All information about a cluster is stored at the 
root site. 

It is technically simple to identify the types of all clusters touching the border: The 
algorithm scans along a border like the one shown in figure 8 and assigns a flag to the 
root of each (occupied) site it is visiting, according to the location of the site visited. 
The flags gathered at each root then identify the borders that the represented clusters 
touch. The P7v(T, n, r) are based on the statistics of these flags. 
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Figure 9. The HK representation of the border of the configuration shown in figure 1 
with the hatched square "activated" to allow for a wrapping cluster. The tuples 
indicate the distance of a site to the site it is pointing to. These numbers are only 
given for the two paths which are relevant for the emergence of a wrapping cluster 
when the system is "glued" along the W and E borders. The dashed line in the top 
right-hand corner indicates the special treatment of a corner, while the dashed arrows 
indicate auxiliary pointers. The highlighted sites give rise to a wrapping cluster. 



Identifying spanning and wrapping 

Starting from a configuration with open boundaries (for example figure 8), spanning 
clusters on a cylinder can be detected easily, by "gluing" the appropriate borders: First 
the roots are identified for each pair of sites, which become nearest neighbours due to the 
new boundary conditions. Preferably, the root of the smaller cluster is then redirected 
to the other, i.e. its label is overwritten by the index of the other root, which, in turn, 
inherits all properties of the overwritten root, such as cluster size or border flags. The 
new pointer resulting from a vertical gluing is shown as a dashed line in figure 8; clusters 
of type U{7v,s} are then spanning on the cylinder. 

It is significantly more complicated to detect wrapping clusters, because these 
clusters cannot be defined by properties of individual sites alone. 

The simplest solution, which is applicable to any kind of topology, such as a cylinder, 
torus or Mobius strip, is to assign to each site an additional set of numbers, which 
indicates its distance to the site it is pointing to. This distance can be measured in 
units of lattice spacings or, even simpler, in multiples of vr: each type of border is 
mapped to a set of integers like 

AT -^(0,1) 5^(0,-1) E^(1,0) W^^(-1,0) (30) 

and for each site an additional flag indicates its location. The distance between two 
sites is then the difference between the tuples associated with their location flag. For 
example, the distance between W and S is (1,-1). In the open system, each site 
is assigned a tupel indicating the distance to the site it is pointing to. The periodic 
boundary conditions are then applied by gluing along the vertical borders. Two clusters 
merge in the way described above, with the additional assignment of a distance vector to 
the redirected root site. This distance vector indicates the distance from the redirected 
root site to the new root site, given by the difference between the distances of the glued 
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sites and their roots. In that way, for any site in a cluster the distance to the root is 
given by the sum over all distance vectors along the tree to the root. 

If the gluing procedure now comes across two sites on either border, which belong 
to the same cluster§, the pathlength to the root is calculated for each of them. If their 
difference is non-zero, a wrapping cluster has been found, the number indicating the 
winding number as a multiple of vr. In the example in figure 9, the left hand path has 
length (1, —1) + (0, 2) = (1, 1), the right hand path (—1, 1), differing by the expected 
length (0,2). It is a topological fact that higher winding numbers cannot appear on a 
cylinder. 

Corners require special attention, because they belong to two different borders at 
the same time. In the example shown in figure 9, all paths shown on the right-hand side 
have length (0, 0), apart from the corner site, which is thought as carrying an "internal 
pointer" to itself, connecting the E border to the border. The internal pointer can 
be understood as follows: if a site points to a corner, it has a connection to two borders 
at once. However, the distance vector of the pointing site can only indicate the distance 
to one border. A convention is required to lift this degeneracy, for example that all 
pointers pointing towards and away from a corner site do so with respect to the or 
S boundary, never with respect to the E oi W boundary. However, it might happen 
that another corner is glued to the corner site, or a neighbouring site is connected to it 
as aW or E site. In this case, auxiliary pointers are introduced: one pointer from the 
connected site to the corner site's W oi E part and another pointer internally connecting 
the corner site's W oi E part to the N or S part. In figure 9 these auxiliary pointers 
and the sub-partitioning of a corner site are shown with dashed lines. 
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